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Abstract 

J> ' We calculate the leading direct-instanton contributions to the operator prod- 

00 ' 

in 

^^ \ up improved QCD sum rules for the nucleon magnetic moments. Remarkably, 

the instanton contributions are found to affect only those sum rules which had 

Q_i| previously been considered unstable. The new sum rules show good stabil- 

ity and reproduce the experimental values of the nucleon magnetic moments 
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uct expansion of the nucleon correlator in a magnetic background field and set 



with values of x, the quark condensate magnetic susceptibility, consistent with 



rN ' other estimates. 
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For over two decades QCD instantons have been associated with fundamental aspects of 
strong interaction physics, as, for example, with the #- vacuum ]TJ, the issue of strong CP 
violation [Q], and the anomalously large rf mass ||. Unequivocal and quantitative evidence 
for their role in hadron structure, however, turned out to be much harder to establish. This 
is mainly due to the complexity involved in dealing with interacting instanton ensembles 
and their coupling to other vacuum fields over large distances. 

Instanton vacuum models || attack the first part of this problem directly, by approx- 
imating the field content of the vacuum as a superposition of solely instantons and anti- 
instantons. This approach has been developed for more than a decade and can describe an 
impressive amount of hadron phenomenology ||. More recently, QCD lattice simulations 
began to complement such vacuum models by isolating instantons in equilibrated lattice 
configurations and by studying their size distribution and their impact on hadron correla- 
tors [§J. While the results obtained with different, currently developed lattice techniques 
have not yet reached quantitative agreement, they do confirm the overall importance of 
instantons and some bulk properties of their distribution in the vacuum. 

Another approach towards linking the instanton component of the vacuum to hadron 
properties has been developed over the last years by generalizing the nonperturbative op- 
erator product expansion (OPE) and QCD sum rule techniques |||-0. While its range of 
applicability is more limited than that of instanton vacuum models and of lattice calcu- 
lations, it avoids the need for large-scale computer simulations and takes, in contrast to 
instanton models, all long-wavelength vacuum fields and also perturbative fluctuations into 
account. Furthermore, the approach is largely model-independent and allows the study of 
instanton effects in a fully analytic and, therefore, rather transparent fashion. 

In the present paper, we adapt this approach to properties which characterize the re- 
sponse of the hadronic system to a weak external field. Specifically, we calculate the direct- 
instanton contributions to the nucleon correlator in the presence of a constant electromag- 
netic field. With the help of background-field sum-rule techniques due to Ioffe and Smilga 
H and Balitsky and Yung ||, we then study previously neglected instanton effects in the 



QCD sum rules for the magnetic moments of the nucleon. 

As already mentioned, our work is based on the nucleon correlation function 

^|rf 4 xe^(0|Tr / (x)r ? (0)|0) F = U (p) + y^^U^(p)F^, (1) 

in the background of a constant electromagnetic field F^. The interpolating fields r](x) with 



proton or neutron quantum numbers |T0|] are composite operators of massless up and down 
quark fields: 



r) p (x) = [u a (x)C la u b (x)] l5l a d c (x)e abc , Vn = Vp (u <-> d). (2) 

For the application in QCD sum rules we need a theoretical description of the correlator ([[]) 
at momenta s = —p 2 ~ lGeV 2 , i.e. at distances x < 0.2 fm. 

The information on the magnetic moments is contained in the second term on the right- 
hand side of Eq. ([!]). It characterizes the linear response of the nucleon to the external field 
and can be decomposed into three independent Lorentz and spinor structures: 

iV(p) = {$0^ + a^f) riiQo 2 ) + i^tyPv - ^vVv)i> n 2 (p 2 ) + o> n 3 (p 2 ). (3) 

Note that the invariant amplitude Hi corresponds to the chirally-even part of the correlator, 
while II2 and II3 are associated with the chirally-odd part. 

The nonperturbative OPE |~T],[8] of the above correlator at small distances can be gener- 



ated by splitting each diagram contributing to (||) in all possible ways into a hard and a soft 
subgraph. The hard subgraphs contribute to the Wilson coefficients and are usually calcu- 
lated perturbatively, with the integration range of each internal momentum restricted^] to be 
larger than the OPE scale fi ~ 0.5 GeV. The soft subgraphs correspond to hadron-channel 
independent condensates, renormalized at \x. In the presence of an external electromagnetic 
field the OPE (up to eight-dimensional operators) involves the additional, Lorentz-covariant 
condensates 



1 In practice, this restriction is often unnecessary (see below). 



(0|g<V9|0> = v^xF»v(0\qq\0), (4) 

g(0\qG^q\0) = v / 4^KF^(0|gg|0), (5) 



^(0|g7sG^g|0) = -V4xaSF^,(0\qq\0). (6) 

(GW = |A a G^, G^ = \e^ up(7 G p(7 with e i23 = -1-) The parameters x, «, and £ play the 
role of generalized susceptibilities and quantify the vacuum response to weak electromagnetic 
fields. The magnetic susceptibility of the quark condensate, \i f° r example, originates from 
the induced spin alignment of quark-antiquark pairs in the vacuum. Note also that \ is 
associated with the lowest-dimensional induced condensate, which enhances its role in the 
OPE and the corresponding sum rules. 

The OPE of ILjj,(p) up to operators of dimension eight, with perturbatively calculated 
Wilson coefficients, has been obtained in Ref. ||. An inherent assumption of this calcu- 
lation - and of the QCD sum rule program in general - is that the short-distance physics 
associated with fields of wavelength smaller than /i _1 is predominantly perturbative. It is 
well known, however, that also strong nonperturbative fields of rather small size exist in the 
QCD vacuum. Instantons, i.e. the finite-action solutions of the classical, Euclidean Yang- 



Mills equation |L2] , are paradigmatic examples of such fields and have a crucial impact on 
the vacuum structure. 

The nonperturbative contributions from short- wavelength fluctuations of quarks and 
gluons around instantons are thus neglected in the standard treatment of the OPE. Their 
relative importance, and hence the justification for approximately disregarding them, de- 
pends both on the instanton size distribution in the vacuum and on the quantum numbers 
of the hadronic channel under consideration. Instantons of smaller (average) size p are ac- 
companied by fluctuations of smaller wavelength, and those contribute more strongly to the 
Wilson coefficients. The hadron-channel dependence of the instanton contributions orig- 
inates mainly from the chirality and spin-color coupling of the quark zero-modes in the 
instanton background. 

Instanton-induced effects are particularly large in the pseudoscalar-isovector and scalar- 



isoscalar channels, because there the quark zero-modes contribute with maximal strength. 
As a consequence, instanton contributions dominate already at short distances in the pseu- 
doscalar sum rules, and are essential for their stability |J. In the vector and axial- vector 
channels, on the other hand, zero-mode contributions are, to leading order in the instanton 
density, absent. 

The strength of instanton contributions to the nucleon channel lies about halfway in- 
between these extreme cases. In the chirally-odd amplitudes, it was found to be roughly of 
the same magnitude as that of the condensate contributions H,0], since the spin-0 diquark 
operators in the interpolating fields (0) couple strongly to instantons. The nucleon channel 
is therefore well suited for studying the interplay between instantons and other vacuum fields 

i- 

The calculation of the leading direct-instanton contributions to the background-field 
correlator ([!]) proceeds essentially along the lines described in Refs. 0-0, to which we 
refer for more details. Similar to the condensates, the bulk properties of the instanton size 
distribution are generated by long-distance vacuum dynamics and have thus to be taken as 
input for this calculation. As before |§-0], we will use the standard values [T3|] p ~ | f m 
for the average instanton size and R ~ lfm for the average separation between neighboring 
(anti)instantons. The results of instanton vacuum models PlJlij confirm these scales, while 



those from the lattice 0] are not yet fully consistent but lie in the same range (with maximal 
deviations of about 50 %). 

Since the average instanton size is of the order of the inverse OPE scale, p < yr 1 ~ 
0.4 fm, instanton corrections to the Wilson coefficients can be substantial. Moreover, since 
p _1 3> Aqcd, these corrections are essentially semiclassical, and at the relevant distances x ;$ 
0.2 fm <C R multi-instanton correlations should be negligible. The instanton contributions 
to the OPE coefficients can therefore be calculated in semiclassical approximation, i.e. by 
evaluating the correlator ([!]) in the background of the instanton and anti-instanton field 
and by then averaging the instanton parameters over their vacuum distributions. (The 
distribution of the instanton's position and color orientation is uniform, due to translational 



and gauge invariance). 

We treat the zero-mode sector of the quark propagator in the instanton field exactly 
and approximate the continuum modes, as before ||, by plane waves. The recently found 
zero-mode dominance of the ground-state contributions to the pion and p-meson correlators 
on the lattice [JTJJ supports the validity of this approximation. The impact of the remaining 
vacuum fields (including other instantons) on the instanton contributions is accounted for 
in a mean-field sense [|l6j and generates an effective mass m(p) = — \^ 2 p 2 {qq) for the quark 
zero modes. We further use the approximate instanton size distribution n(p) = n5(p — p) fy 



and the self-consistency condition JT7| (qq) = —2 J dp ^4 = —2 -^y=y (which is numerically 



satisfied to good accuracy) to eliminate the n dependence from the resulting expressions. 

Up to operators of dimension eight, we find the leading instanton contributions to the 
correlator ([!]) to arise from just one type of graph, in which two of the quarks (emitted 
from the current (|2|) at x = 0) propagate in zero-modes while the third interacts with the 
background field through the magnetized quark condensate. Graphs in which the back- 
ground field couples directly to a hard quark in a zero-mode, vanish. The same holds for 
graphs in which the interaction with the photon causes a transition from zero-mode to 
continuum-mode propagation?]. The contribution from direct instantons is thus generated 
by the interplay between the rather localized quark zero modes and more slowly varying, 
nonperturbative vacuum fields. Effects of this type appear naturally in the OPE, whereas 
they are difficult to account for in, e.g., quark models [0. 

Evaluating the corresponding graphs in the instanton and magnetic background fields as 
described above, and averaging over the instanton size distribution, we arrive at 

(01^(^^(0)10)^ = nr(z) - ^^(^q) F cv / d 4 x (r2 - _ 2) l {xl - _ 2)3 (7) 

(r = x — Xq, where Xq specifies the center of the instanton) for the proton. The corresponding 



2 Contributions of this type are essential in the pseudoscalar three-point correlator associated with 
the pion electromagnetic form factor H. 



neutron correlator is obtained by replacing e u with e^. In order to put the instanton contri- 
bution (0) to use in the sum rules, we also need its Fourier and Borel transform^. Again for 
the proton, it reads 

n 3 (M 2 ) = -5L- aX ?M*I(*>). (8) 

Here, M denotes the Borel mass parameter. We have also used the standard definitions 
z = Mp, a = —(2ir) 2 (qq), and abbreviated the integral 



1 ( z ) = / ~TTi v^ e 4Q(1_Q) = 4e 2 

jo cr(l — a) z 



z 2 \ „ (z T 



Mt +K ^2 



(9) 



Note that (P) shows the typical exponential Borel-mass dependence of instanton contribu- 
tions [||]. Together with the appearance of the new scale p, this distinguishes them from the 
logarithms and power terms of the standard OPE. 

An important qualitative property of the instanton contribution has been made explicit 
in Eq. (|j): to leading order, direct instantons contribute almost exclusively^ to one invariant 
amplitude, II3, which is associated with the chirally-odd Dirac structure a^ v . Exactly this 
amplitude was singled out in the previous sum rule analysis of Ref. 0, for two reasons: first, 
one of its Wilson coefficients contains, in the "pragmatic" version of the OPE (see below), 
an infrared divergence. Secondly, the II3 sum rule failed to show a fiducial Borel region of 
stability Jl8| (even after proper subtraction of the divergence), while the other two sum rules 
provide stable and accurate values for the nucleon magnetic moments even without direct- 
instanton corrections, and with a value of \ consistent with other, independent estimates 

El- 



3 For the Borel transform we follow the convention of Ref. JR 

4 There is also a small direct-instanton contribution of similar structure to the amplitude n2. 
This term turns out to be too small to have an appreciable impact on the corresponding sum rules, 
however, and will not be discussed further. 



To understand the first point we note that, since in QCD perturbative contributions 
from soft loop momenta are normally small compared to the corresponding condensate con- 



tributions, it is standard procedure not to remove them in sum-rule calculations |]20|| . This 
simplification - which goes under the name of "pragmatic OPE" - fails, however, if infrared 
divergences appear in diagrams associated with a Wilson coefficient. Such an infrared di- 
vergence was encountered in the OPE of II3, in a graph where a vacuum gluon field and the 
background photon interact with the same hard quark line. Hence the contribution from 
soft loop momenta has to be cut off explicitly, according to the rules of the exact OPE. A 
similar divergence was found before |2T| in the vector meson correlator when two soft gluon 
fields couple to the same quark line. 

The authors of || conjectured that the appearance of infrared singularities in II3 and the 
absence of a stability region in the associated sum rule might be related. This seems unlikely, 
however, in view of the later finding [^] of similar infrared divergencies (which always occur 



when a quark line interacts with multiple soft gauge fields) in the other two amplitudes of 
(0), which nevertheless lead to satisfactory sum rules. Direct instantons, on the other hand, 
contribute almost exclusively to II3, and their previous neglect could offer a more plausible 
explanation for the instability of the II3 sum rule. Support for this conjecture, which we are 
going to test quantitatively below, comes from two other chirally-odd nucleon sum rules (for 
the nucleon mass || and its isospin splitting |7|), where exactly such a selective stabilization 
due to instantons has been found. 

Our modified II3 sum rule is obtained by equating the Borel transform of the standard 
OPE from Ref. || and the instanton contributions (H) to the Borel-transformed double 
dispersion relation for the correlator (P, with a spectral function parametrized in terms of 
the nucleon pole contribution (containing the magnetic moments) and a continuum based on 
local duality. Including the infrared- divergent term encountered in Ref. M, duly truncated 



at the OPE renormalization point, the new sum rule for the proton readsQ 
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(10) 



where m is the nucleon mass, W the continuum threshold, and X^ the coupling of the 
current @ to the nucleon state, (0\r]\N) = X^u. For the mixed quark condensate we use 
the standard parametrization (qa^G^q) = —m^(qq) with m 2 , = 0.8 GeV 2 , and Jem — 0.577 
is the Euler-Mascheroni constant. The additional parameters A pn determine the strength 
of electromagnetically induced transitions between the nucleon and its excited states. The 
sum rule for the neutron is obtained from (|10"D by interchanging e u and e d and by replacing 
Up, /i" — >■ n n and A p — >■ A n . We have also defined A^ = 327r 4 A 2 v , L = ln(M/A)/ln(/x/A) 



(A = 0.1 GeV), and transferred, using the standard expressions 

i fw^y 



E n (M) = l-e 



l + \ ' j\ \M 2 J 



(11) 



the continuum contributions to the OPE-side of the sum rules. The appropriate form of 
these contributions has recently been clarified in Ref. [^ . 

In principle, several alternative options are available for the quantitative analysis of 
background-field sum rules. In practice, however, one is limited by the fact that their fiducial 
domain (i.e. the Borel-mass region in which the neglect of higher-order terms in the short- 
distance expansion is justified while the nucleon pole still dominates over the continuum) 
is generally not large enough to determine all the unknown parameters from a stable fit. 
The authors of Ref. succeeded, however, in eliminating the susceptibilities and other 
constants by combining the two sum rules for Hi and n 2 and their M 2 -derivatives. The 
magnetic moments can then be fitted and are in good agreement with experiment, although 
taking derivatives of the sum rules generally reduces their reliability. Unfortunately, the 



'We have corrected some errors which appeared in the expressions of Ref. S 



above procedure also eliminates the continuum contributions, which makes it impossible to 
determine whether a fiducial stability domain exists. 

In any case, this procedure is ineffective for our new sum rules (|10|), as it cannot eliminate 
the additional ^-dependence introduced by the direct-instanton contribution. Similarly, 
working with the ratio of background-field and mass sum rules, as done in Ref. J23J for 
the 111 sum rules, offers no advantage in our case since already the leading terms in the 
OPE of the chirally-odd mass and background-field sum rules differ. Therefore, we resort 
to a direct minimization of the relative deviations between the two sides of Eq. fllPf) . The 
coupling A 2 = 2.93 GeV 6 and the continuum threshold W = 1.66 GeV are, following the 
procedure of Ref. || , obtained by fitting the instanton- improved nucleon mass sum rule || 
to the experimental nucleon massQ. The values of the two susceptibilities k = —0.34 ±0.1 



and £ = —0.74 ± 0.2 were estimated in independent work by Kogan and Wyler [25]. This 
enables us to fit both sides of the sum rules ( |T0D by varying x an d A p (or A n , respectively) 
while keeping the magnetic moments fixed at their experimental values. 

The fits are performed in the fiducial Borel mass domain, which is bounded from below 
by requiring the highest-dimensional operators to contribute at most 10 % to the OPE and 
from above by restricting the continuum contribution to maximally 50 %. The resulting 
fiducial domains of both the proton and neutron sum rules, 0.8 GeV < M < 1.15 GeV, are 
larger than those of the sum rules based on Hi and n 2 . 

The fit results are shown in Fig. [I], where for both the proton and the neutron sum 
rules the direct-instanton contributions, the remaining OPE including the continuum con- 
tributions, their sum (which makes up the left-hand side of Eq. (|10D), and the right-hand 
sides are plotted. The fit quality of both the proton and neutron sum rules is excellent. As 
previously in the instanton- improved, chirally-odd nucleon mass sum rule pfl, the theoretical 



6 We have included anomalous-dimension corrections in the mass sum rule of Ref. || and restored 
the full four-quark condensates, since we are dealing exclusively with Ioffe's interpolating field (|2|). 



10 



side of the sum rules, including the instanton- induced part, is almost indistinguishable from 
the phenomenological side. Moreover, Fig. [I] shows that the direct-instanton contributions 
can reach about half the magnitude of the remaining terms in the OPE, which makes it 
evident why their previous neglect had a detrimental impact on the stability properties. 

An alternative way of evaluating the optimized sum rules consists in solving them for /xjv 
and plotting the result as a function of the Borel mass, as shown in Fig. ^|. The resulting 
functions /x Pin (M) specify the value of the magnetic moment which is required to make both 
sides of the sum rule fllQl) match exactly at each value of the Borel mass. The instanton- 
corrected sum rules render fi(M) practically M-independent, thereby again indicating an 
almost perfect fit for both proton and neutron. 

This fit predicts x — —4.96 GeV -2 for the proton and x — —4.73 GeV -2 for the neutron 
sum rule. These values correspond to the OPE scale \x = 0.5 GeV adopted for our sum 
rules and lie inside the range obtained from other estimates f8,|l9")22|1. They are somewhat 



smaller in magnitude than the value x — —5.7 GeV 2 found in the two- and three-pole 
models of Ref. [BJ|. (Our predicted values for the excited-state transition parameters are 



A p ~ 0.28GeV 2 and A n ~ -0.27GeV~ 2 .) 

In conclusion, we have recovered a third reliable sum rule for the nucleon magnetic 
moments. In contrast to the other two, it receives previously neglected direct-instanton 
contributions which arise from the interplay with long-wavelength vacuum fields. Our new 
sum rule is built on the chirally-odd amplitude n 3 of the nucleon correlator in an electro- 
magnetic background field and found to be at least as stable as the other two, although it 
had previously been regarded as flawed. The new sum rule adds to the predictive power of 
the background-field sum rules and strengthens their mutual consistency. 

Furthermore, our results reinforce a systematic pattern which emerged from previous 
studies of direct-instanton effects both in the nucleon and pion channels: those sum rules 
which worked satisfactorily without instanton corrections receive little or no direct instanton 
contributions, and previously less reliable or completely unstable sum rules are stabilized 
by large instanton contributions. This pattern points not only towards the importance of 

11 



direct instantons in particular sum rules, but also supports the adequacy of their semiclassical 
implementation into the OPE. Our results show that these conclusions continue to hold in 
the presence of a "magnetized" vacuum. 

M.A. and M.K.B. acknowledge support from the U.S. Dept. of Energy under grant 
number DE-FG02-93ER-40762. H.F. acknowledges support from Deutsche Forschungsge- 
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FIGURES 
FIG. 1. The OPE (dashed line) and direct instanton (dotted line) contributions to the new u^v 

sum rules for the proton (positive range) and neutron. Their sum (dot-dashed line) is compared 

to the RHS (solid line). 

FIG. 2. The Borel mass dependence of the magnetic moments of the proton (upper) and neutron 
calculated from the optimal fit of LHS and RHS. 
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